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ABSTRACT 


Raman-based spectroscopy is potentially militarily useful for standoff detection of 
high explosives. Normal (non-resonance) and resonance Raman spectroscopies are both 
light scattering techniques that use a laser to measure the vibrational spectrum of a 
sample. In resonance Raman, the laser is tuned to match the wavelength of a strong 
electronic absorbance in the molecule of interest, whereas, in normal Raman the laser is 
not tuned to any strong electronic absorbance bands. The selection of appropriate 
excitation wavelengths in resonance Raman can result in a dramatic increase in the 
Raman scattering efficiency of select band(s) associated with the electronic transition. 
Other than the excitation wavelength, resonance Raman is performed experimentally the 
same as normal Raman. In these studies, normal and resonance Raman spectral 
signatures of select solid high explosive (HE) samples and explosive precursors were 
collected at 785 nm, 244 nm and 229 nm. Solutions of PETN, TNT, and explosive 
precursors (DNT & PNT) in acetonitrile solvent as an internal Raman standard were 
quantitatively evaluated using ultraviolet resonance Raman (UVRR) microscopy and 
normal Raman spectroscopy as a function of power and select excitation wavelengths. 
Use of an internal standard allowed resonance enhancements to be estimated at 229 nm 
and 244 nm. Investigations demonstrated that UVRR provided '-'2000-fold enhancement 
at 244 nm and -SOO-fold improvement at 229 nm while PETN showed a maximum of 
-'25-fold at 244 nm and -'190-fold enhancement at 229 nm solely from resonance effects 
when compared to normal Raman measurements. In addition to the observed resonance 
enhancements, additional Raman signal enhancements are obtained with ultraviolet 
excitation (i.e., Raman scattering scales as for measurements based on scattered 
photons). A model, based partly on the resonance Raman enhancement results for HE 
solutions, is presented for estimating Raman enhancements for solid HE samples. 
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I. 


INTRODUCTION 


A. MOTIVATION 

The development of reliable and effective explosive detection technologies is an 
important challenge for military and homeland security applications. From the wide 
proliferation of landmines, termed the poor man’s artillery, to the pressing need to detect 
improvised explosive devices, the need for a dependable and effective explosive 
detection system is becoming increasingly important. Due to the extreme challenges that 
hinder effective explosive detection, it is widely recognized that there is no “silver bullef ’ 
solution. Several challenges hinder the development and deployment of standoff 
explosive detection systems. Landmines and improvised explosive devices (lED) are 
intentionally concealed explosives designed to thwart detection efforts. Additionally, 
most commercial and military grade explosives exhibit extremely low vapor pressures, 
which significantly limits the ability to detect explosives in the gas phase at standoff 
distances or otherwise. Every explosive detection system developed or theorized to date 
has drawbacks when evaluated against the following factors: minimum detectable limit of 
explosive concentrations, probability of detection, probability of false alarms, the variety 
of explosives that can be reliably detected, portability of the unit, analysis time, standoff 
capability, container or soil penetration capability, required operator training, and the cost 
of the system. 

Due to the numerous and ever changing number of ways that explosives are used 
against military personnel and equipment, it is important to develop a suite of explosives 
detection technologies. Eor example, technologies developed to detect landmines may be 
unsuitable to detect explosives on a suicide bomber. The chemistry and physics of 
explosive detection leads to cases where one technology is very effective within inches of 
the explosives (e.g., used in personnel screening) but would be completely ineffective at 
standoff distances. For these reasons a layered and multifaceted approach to landmines 
and lED defense is necessary. 
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Currently, there are a variety of teehniques used to deteet explosives: sniffing 
explosive vapors (e.g., fluoreseenee quenehing and ion mobility speetroseopy), imaging 
bulk explosives (e.g., x-ray baekseattering and millimeter wave imaging), and 
speetroseopie deteetion of explosives (e.g., infrared, Raman, and laser-indueed 
breakdown speetroseopy teehniques, and nuelear teehniques). Due to the drawbaeks 
experieneed by all explosive deteetion systems, it is highly desirable to eombine 
deteetion systems in orthogonal approaehes [1], meaning deteetion methods should 
signifieantly differ in deteetion modality to eapitalize upon the strengths of individual 
teehniques while minimizing overall system weaknesses. When orthogonal systems are 
used optimally it is possible to obtain a probability of deteetion and false alarm rate that 
is more favorable than any individual deteetion modality. 

This researeh foeuses on fundamental aspeets of resonanee Raman speetroseopy 
as applied to seleet HE materials to advanee our eurrent understanding and potentially 
improve HE deteetion at standoff distanees. 

B. BACKGROUND OF PREVIOUS WORK 

Previous researeh into standoff deteetion of explosives using Raman speetroseopy 
eondueted by J. C. Carter’s researeh group at Eawrenee Eivermore National Eaboratory 
and S. K. Sharma’s researeh group at the University of Hawaii provide the foundation of 
this investigation. Eaeh of the aforementioned researeh groups has eompleted 
experiments that demonstrate the initial feasibility of using normal Raman speetroseopy 
to deteet explosives at standoff distanees; however, it is a mueh greater ehallenge to 
demonstrate a teehnology’s military utility sinee this involves proposed eoneepts of 
operations, theorizing system required performanee parameters, and utilizing more 
representative target sets. 

The work eompleted in this thesis exploits the previous results derived from work 
from several researeh groups investigating standoff Raman and resonanee Raman to 
further improve the ability to deteet high explosives at standoff distanees by assessing the 
potential gains aehievable using resonanee Raman to deteet explosives. 
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C. RESEARCH GOALS 

The primary goal of this research is to demonstrate and quantify the 
improvements obtainable using resonance Raman spectroscopy to detect high explosives. 

D. RESEARCH OBJECTIVES 

The primary research objectives were developed to achieve the stated research 
goal. The key objectives include the following; (1) performing absorption studies of 
select HE materials (e.g., TNT, PETN, PNT, DNT) to determine absorption bands 
throughout the UV-Vis spectral region for acquiring UV Raman spectral signatures; (2) 
performing solution-based, UV-excited, Raman HE measurements while utilizing an 
internal Raman standard measured in-situ; (3) determining the resonance Raman signal 
enhancements, if any, from the solution-based HE measurements relative to normal 
Raman (4) investigate experimental parameters in terms of two different UV 
wavelengths, and select laser power settings. 
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II. DETECTION OE EXPLOSIVES 


A. SELECT PROPERTIES OF EXPLOSIVES 

Explosives are unstable moleeules that ean share many similar physieal 
eharaeteristies. Most military explosives tend to have a melting point greater than 80°C, 
have densities typieally between 1.6 and 2.0 g em'^, and they are thermally stable and 
relatively insensitive to shook and other stimuli [2]. High explosives have low densities 
sinoe they are oomprised of low atomio number (Z) atoms like oarbon, nitrogen oxygen, 
and hydrogen. High explosives, with some notable exoeptions suoh as peroxide-based 
explosives, oontain one or more nitro funotional groups (-NO 2 ) that tend to be the least 
stable bond in the explosive moleoule [3]. While most explosive oompounds oontain 
nitro groups, the organio ohemistry of explosives ean vary widely. Explosive moleeules 
ean be olassified into several organio ohemistry groups inoluding open and olosed ohain 
nitro-aliphatio, nitrate esters (e.g., NC, NG, PETN), and nitramines (e.g., RDX, HMX) 
and aromatio oompounds (e.g., DNT, TATB, TNT). 

I. Types of Explosives and Chemical Composition 

Due to stability and availability, oommeroial and military grade explosives are the 
types of explosives most frequently enoountered on today’s battlefield. The elemental 
oomposition of oxygen and nitrogen is unusually high in most military explosives 
oompared to other innoouous oompounds, an exploitable feature for some analytioal 
deteotion teohniques. The homemade explosive triaoetone triperoxide (TATP) is at the 
top of the list of notable exoeptions to this rule of thumb; while it oontains a high 
peroentage of oxygen it does not oontain nitrogen atoms. The average nitrogen and 
oxygen oontents of explosives are 31 ± 12% nitrogen and 45 ± 8% oxygen [4]. 
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Figure 1. Weight percent of oxygen and nitrogen in select high explosives (From; [4]). 


2. Peroxide-based Explosives 

In 2001, the U.S. encountered homemade peroxide-based explosives when 
Richard Reid attempted to detonate PETN (pentaerythritol tetranitrate) with TATP on a 
Miami bound Boeing 767 flight. TATP was also used in the London public transit 
system bombings in 2005, and it has been a popular choice among homicide bombers for 
decades since it is very easy to synthesize. Before Operation IRAQI FREEDOM, TATP 
was one of the most common constituents used in lEDs, since it is effective and could be 
manufactured with easily obtained ingredients. The major drawback with TATP is its 
instability, which makes handling and synthesis dangerous. When factors allow groups 
to readily obtain commercial or military explosives, these will be preferred over 
peroxide-based explosives for performance and safety reasons. TATP has a significantly 
higher vapor pressure than most military explosives; for a given mass, TATP has 
approximately 10^ more molecules in the air at room temperature than TNT (2,4,6- 
trinitrotoluene or also 2,4,6-TNT) [4]. 
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B. SPECTROSCOPY OF EXPLOSIVES 


Spectroscopy was traditionally defined as the study of interaction of 
electromagnetic radiation with matter as a function of wavelength [5]. This includes the 
transmission, absorption, scattering and emission of light from the sample. This 
definition was later expanded to encompass other forms of analytical measurements that 
include particles interacting with matter. After the analyte is perturbed, several different 
observable phenomena may be detected to provide characteristic information concerning 
the probed atom or molecule. Spectroscopic measurements may be based on absorption, 
scattering, or characteristic photon emissions. Several different forms of spectroscopy 
have been developed to probe solid, liquid and vapor phase explosives to provide 
qualitative and quantitative information. A recent publication by D. S. Moore [6] 
provides a thorough review on this topic. 

1. Explosive Vapors 

Most high explosives have exceedingly low vapor pressures (see Figure 2), so 
very few molecules in the gas phase are present outside the solid or liquid explosive 
device. The concentration of most military grade high explosives present in the gas 
phase at room temperature can be measured in parts per trillion (ppt) to parts per million 
(ppm). As an example, converting this concentration to more familiar units yields 
approximately 0.1 ng cm' for 2,4,6-TNT [7]. Some homemade peroxide-based 
explosives, like triacetone triperoxide (TATP), have considerably higher vapor pressures 
that can be measured in parts per thousand. 
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Figure 2. Vapor Phase Coneentration of eornmon explosives as a funetion of 
temperature. The solid lines indieate experimental values (From; [6]). 

In addition to the extremely low eoneentration of explosive moleeules present in 
the gas phase under normal eonditions, nitrogen-based high explosives have a tendeney 
to rapidly eondense onto aerosols or solid surfaees rather than remain in the vapor phase 
further hindering deteetion efforts. 

Canines have an exeeptionally aeute sense of smell and have typieally been the 
gold standard for explosive vapor deteetion, however some traee deteetion systems have 
reeently demonstrated improved sensitivities over eanine deteetion by deteeting sub- 
pieogram quantities of seleet explosives [8]. Even with this exeeptional sensitivity, it is 
diffieult to deteet explosive vapors using a point sensor at distanees greater than a few 
meters. Sometimes favorable wind eonditions ean extend the range of explosive vapor 
deteetion, however, loealizing sueh sourees from turbulent plumes over extended 
distanees with point sensors has not been demonstrated to date. 

The presenee of explosive related materials (ERM) sueh as tagging agents, 

impurities, reaetion byproduets, plastieizers or deeomposition produets or enhaneers sueh 

as propane in explosive deviees may be more traetable than direet deteetion of the gas 
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phase explosive moleeules sinee ERM and enhaneers tend to have mueh higher vapor 
pressures. The United Nation’s International Civil Aviation Organization (ICAO) 
requires the addition of tagging agents to high explosives to aid in explosive deteetion 
[9]. Tagging agents are added in quantities varying from 0.1 to 0.5% by weight. 2- 
nitrotoluene and 4-nitrotoluene (or p-nitrotoluene or PNT) are two of the four eommonly 
used taggants to faeilitate improved vapor deteetion of explosives. Natural impurities left 
over from ehemieal synthesis reaetions typieally aid explosives deteetion efforts. For 
example, 2,4-dinitrotoluene (2,4-DNT) is found in TNT due to ineomplete nitration 
reaetions. DNT has a higher vapor pressure than TNT and may be more easily deteeted 
than TNT itself [1,9]. 

2. Bulk Explosives 

A few speetroseopie methods have been developed to probe bulk explosives. 
Sinee explosives in landmines or lEDs are eoneealed by soil or in vehieles or under 
elothing, bulk speetroseopy methods must be eapable of penetrating multiple inehes of 
ooneealing material. Direet line-of-site deteetion teehniques, sueh as Raman and EIBS 
(laser-indueed breakdown speetroseopy), are very limited in sueh seenarios beeause 
signal transduetion into and out of the deviee is often limited. Bulk deteetion methods 
typieally have limited sensitivities that make it diffieult or impossible to deteet small 
antipersonnel mines that eontain 50 to 100 g of high explosives, and the exeitation souree 
typieally must be less than a few meters from the target. Even so, there are a few bulk 
deteetion methods sueh as pulsed fast and thermal neutron analysis (PFTNA) that are 
useful for bulk explosives deteetion. PFTNA uses energetie neutrons to probe high 
explosives though several layers of ooneealing materials to oause oharaoteristio gamma 
ray emissions that oan be deteeted by sointillators [10]. In this method, and similar to 
other bulk deteetion methods, the neutron souree must be within inehes or at most a few 
meters from the target to make a qualitative determination of the presenoe of explosives. 
Nuolear quadrupole resonanoe (NQR) is another bulk deteetion method proven effeotive 
in handheld landmine deteetion at olose distanoes by probing the magnetio dipole 
moments of in bulk explosives [11]. 
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3. 


Surface Residues 


The surfaces of conventional and improvised munitions are contaminated with 
explosive residues [12, 13]. While concentrations of these explosive contaminants can 
vary, trace quantities can be readily detected by trace detection techniques (e.g., ion 
mobility spectroscopy and fluorescence quenching). Additionally, the transfer of 
explosive residues while handling explosives is typically unavoidable. Personnel 
involved in the storage, transportation and emplacement of mines and lEDs typically 
have identifiable amounts of explosive residues on their hands, clothes and personal 
effects, and this contamination can extend to vehicles, residences and other items that 
come into contact with the contaminated individual [14]. Explosives are very polar 
molecules and tend to strongly physisorb to many substrates. Most military grade 
explosives can be detected on persons who handle explosives over 48 hours past the time 
of exposure even when vigorous washing with surfactants is performed [13]. The typical 
quantity of explosives found in the fingerprints of contaminated individuals is 
approximately 40 to 100 //g of material [7, 14]. Raman microscopy (table-top device) 
measurements of explosive residues are used in forensic investigations [30]. 

One must consider the preexisting background levels of explosives and 
environmental degradation of explosive residues as important analytical considerations 
when detecting explosive residues. 

4. Other Detectable Characteristics 

While this thesis focuses on spectroscopic detection of explosives using Raman- 
based spectroscopy, there are other notable characteristics of explosives that can be 
exploited. These characteristics include dielectric constants and density. X-ray 
absorption and backscatter techniques have been successfully employed at close ranges 
by exploiting the characteristic densities of high explosive compounds. 

While direct standoff detection of the explosive constituents is considered the 
‘Holy GraiT of landmine and lED detection, other non-direct detection modalities can be 
employed that rely on the device’s metal content, radar signature, placement (e.g.. 
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disturbed soil), seismic signature, electronic signature, infrared heat signature, or 
predetonation. It is commonly recognized that an effective explosive detection system 
will employ orthogonal detection modalities or methods that rely on independent 
measurements [1]. The use of orthogonal approaches can improve the probability of 
detecting explosives while reducing false alarm rates and protecting against simple 
countermeasures. 

C. REQUIRED SYSTEM PERFORMANCE 

Fielding an effective and reliable tactical standoff explosive detection system is a 
daunting challenge. The basic concepts that must be addressed are presented below. 

1. Limit of Detection (LOD) 

The limit of detection (LOD) for a given sensor is a critical figure of merit. In 
spectroscopy, the LOD is typically defined as an analyte signal that has a signal-to-noise 
ratio of two or three depending on the confidence interval of the measurement [1]. The 
ideal tactical standoff explosive detection systems is capable of detecting bulk explosives 
through the concealing materials, detecting extremely low concentrations of explosive 
vapors, or detecting trace explosive residues. Normal and resonance Raman 
spectroscopy lack the ability to detect explosive molecules in the vapor phase at standoff 
distances. Additionally, the direct-line-of-site requirement for Raman measurements is 
severely limiting for explosives detection scenarios requiring the penetration of 
concealing materials. The most likely scenario for standoff Raman-based detection 
involves direct-line-of-site measurements of HE or explosive related material residues. 

2. Time per Analysis 

The time per analysis in Raman spectroscopy is the total required time to 
interrogate the sample via laser excitation, to detect the return analyte signal at a 
specified confidence, and to process the data and report. For forensic or intelligence 
gathering purposes, a time per analysis of a few minutes may be suitable; however, 
tactical combat situations require rapid detection and identification. Ideally a standoff 
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detection method with very rapid time per analysis is employed to allow continuous 
sampling from a moving vehicle at some tactical rate of advance. In cases where the time 
per analysis is more protracted, the detection system may be more suitable as a 
confirmation tool that is employed after visual detection or other explosive detection 
means are employed. 

3. Specificity 

The specificity of an analytical method is the ability of the instrument to correctly 
identify the desired analytes in a representative sample matrix [15]. Ideally, the system 
must be capable of detecting the analyte in the presence of interferences while not 
triggering a false positive alarm. Interferences vary and may include similar chemical 
moieties that produce signatures similar to that of the analyte. For example, an explosive 
detector ideally will not alarm on the presence of nitrogen containing fertilizers or other 
innocuous chemical compounds that may be routinely encountered, otherwise, the 
detector’s effectiveness may be diminished. Convolution of spectral peaks can also occur 
due to matrix effects, which may hinder or preclude explosive detection efforts. As an 
example, standoff Raman-based measurements are capable of detecting O 2 and N 2 in air, 
which could obscure or completely mask the Raman band of a target analyte. 

4. Robustness 

The robustness of an analytical technique is always an important consideration; 
however, in the military this metric is tested to extremes. It is a measure of the ability of 
the analytical method to retain its effectiveness despite intentional perturbations to the 
system [15]. For example, the ideal military explosive detectors should be equally 
capable of detecting explosives in -20°F temperatures, as well as temperatures reaching 
130°F. Additionally, the system should remain effective despite changes in relative 
humidity, precipitation, ambient light conditions, and changes to the sample matrix (i.e., 
the explosive may be located in different matrices and/or on different substrates, etc.). In 
military explosive detection operations, the matrix and/or substrate for a given sample 
will most likely not be known a priori and background subtraction from these and other 
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interferences may not be possible or feasible. Generally speaking, Raman measurements 
of real-world samples are always background noise limited. 

Finally, detection systems must be designed with the warfighter in mind. Such 
systems must conform to the warfighters tactical situations in terms of robustness and 
ease of use including training requirements. Soldiers on the move should not be 
burdened with performing spectral matches even if technical reachback capabilities are 
available. The system should simply be capable of signaling an alarm to the presence of 
explosives at a given confidence level, localize the detection event, and if desired indicate 
the type of explosive that was detected with an indicated confidence level. 
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III. FUNDAMENTALS OF RAMAN SPECTROSCOPY 


A. SCATTERING PROCESSES 

Incident light on a surface will be scattered, absorbed, or pass directly through the 
material. Historically the study of light has inspired much of the theoretical and 
experimental work performed in physics. The nature of light scattering in the sky 
inspired Rayleigh’s important contributions to the nature and behavior of light [16]. Both 
Lord Rayleigh and Sir Raman made fundamental discoveries in physics inquiring about 
the scattering of light. 

1. Raleigh Scattering 

There are several factors that determine the type of scattering that is observed 
when light interacts with matter, but the main factors include the refractive index of a 
particle on a substrate and the size of a particle compared to the wavelength of incident 
light [5]. Rayleigh scattering is an elastic scattering process that occurs when the 
diameter of the particles (the sample have much smaller dimensions than the 

wavelength (T) of the incident electromagnetic waves <0.05/1) [5]. 

For a given Raman measurement, Rayleigh scattering is orders of magnitude more 
intense than the Raman scattered light. For this reason, Rayleigh scattering must be 
filtered before reaching the Raman detector. 

2. Raman Scattering 

The deep blue colors of the Mediterranean Sea observed by Sir Raman in the 
early 1920s inspired spectroscopy discoveries in a newly discovered scattering 
phenomenon now known as the Raman effect [16]. While Rayleigh scattering, centered 
at the incident wavelength, is the dominant scattering effect for the interaction of light 
with matter, a weaker, shifted Raman effect is also observed. The Raman effect is 
significantly weaker than Rayleigh scattering and only involves one Raman scattering for 
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approximately every 10^ to 10'° incident photons [17, 18]. The oscillating electric field 
of the incident laser excitation interacts with matter; a short-lived distortion of a 
molecule’s electron cloud perturbed by vibrations of the molecule gives rise to the 
Raman effect. The Raman effect will bathochromic (red) shift and hypsochromic (blue) 
shift light by two separate processes known as Stokes and anti-Stokes scattering, with 
Stokes scattering being the dominant effect. In Stokes scattering, a molecule in the 
ground vibrational state is excited to a virtual state (i.e., not a real quantum state in the 
molecule), which nearly simultaneously emits a less energetic red-shifted photon. In 
anti-Stokes scattering, a molecule in an excited vibrational state imparts energy to the 
incident photon and relaxes back to the ground vibrational state by emitting a blue-shifted 
photon. The energy level diagram shown in Figure 3 depicts the Rayleigh, Stokes and 
anti-Stokes scattering processes with the m state representing the ground state and the n 
state representing an arbitrary excited state. 


Virtual States 


Vibronic States 



Figure 3. Energy level diagram of normal inelastic Raman processes and Rayleigh 
scattering (After: [18]). At ambient temperature explosive molecules are mostly 
in the ground state (m) and the Stokes effect dominates the Raman scattering 
processes. Rayleigh scattering is an elastic scattering effect. In resonance Raman 
techniques, which will be discussed later, the molecule is excited to an electronic 

transition. 


At ambient temperatures, most molecules are in the ground state and therefore. 


Stokes scattering intensity is much stronger than anti-Stokes scattering [18]. Due to 
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thermal energy, a very small quantity of moleeules will be exeited to higher vibronie 
states at room temperature. Stokes shifted Raman bands are used for analytieal purposes 
sinee these bands are more intense than anti-Stokes bands. 

3. Raman Spectroscopy 

Raman spectroscopy systems are generally comprised of an excitation source 
(e.g., laser), receiver, spectrometer, detector and computer. Laser sources stimulate the 
Raman process, fiber optics, and a microscope or a telescope collects the Raman signals. 
A laser rejection filter is typically used to reject the laser line (i.e., Rayleigh scattered 
light) and a spectrometer disperses light along the focal plane of a detector, such as a 
charge-coupled device (CCD). 


telescope scatter 



Figure 4. Illustration of a fiber-optic coupled standoff Raman detection system (From; 

[19]). The laser beam is coaxial to the telescope field of view. LI is a 
collimating lens and fr is the f-number of the telescope. 

Raman spectroscopy is useful for solids, liquids, and gas phase measurements and 
is currently a widely used analytical technique [18]. There are many analogs between 
Raman spectroscopy and infrared spectroscopy. The unique absorption bands observed 
in infrared and Raman spectroscopy both correspond to molecular vibrations, although 
Raman spectroscopy indirectly probes molecular vibrations. These vibrations can be 
mapped directly to in- and out-of-plane stretches and vibrations of molecular bonds, 
which are excited at characteristic frequencies. Since molecular vibrations are sensitive 
to the local chemical environment surrounding the bond, the infrared and Raman spectra 

forms a unique molecular “fingerprint” that can be exploited. 
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Raman Band 

Frequency (cm'^) 

Band Strength 

CH 3 symmetric stretch 

2862-2882 

Very strong 

C—C stretch 

1040-1100 

Strong 

R-NO 2 asymmetric stretch 

1530- 1600 

medium - weak 

R-NO 2 symmetric stretch 

1310-1397 

very strong 

Aromatic rings 

1600 

very strong 


Table 1. Table of select vibrational frequencies important in the analysis of high 

explosives (After; [17]). 

Raman and infrared spectroscopy measurements are often used to identify 
unknown molecules by comparing the unique spectrum of an unknown to molecular 
databases and vibrational correlation charts. While infrared spectroscopy is an absorption 
phenomenon, the Raman effect is observed by the scattering of light. Raman 
spectroscopy is a weak phenomenon, as previously discussed, and can suffer from 
fluorescence interference due in part to the sample, the sample matrix/substrate, or due to 
laser-induced sample degradation. 

Anti-Stokes lines are much weaker than Stokes lines by a factor of 
exp(-hry / kT ), where h is the reduced Planck constant, co is the angular frequency, k is 
Boltzmann’s constant, and T is temperature in Kelvin. Thus, anti-Stoke lines are 
normally omitted in published Raman spectra unless special anti-Stokes techniques are 
being used. 

Figure 5 shows a sample Raman spectrum illustrating Rayleigh, Stokes, and anti- 
Stokes scattering. The y-axis is the Raman signal intensity and the x-axis is typically 
presented as the Raman shift in wavenumbers (cm'^). 
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Figure 5. Raman spectra of chloroform (CHCI 3 ) obtained using a 514.5 nm excitation 
laser (From; [17]). Note that the Anti-Stokes peaks (left side of the spectrum) are 
much weaker than those obtained from Stokes processes (right side). The strong 
center Rayleigh band is blocked in typical measurements so the weak Raman 
signal is not swamped by the Rayleigh scattered light. 


B. ENHANCING RAMAN MEASUREMENTS 


1. Raman Frequency Dependence 

One of the practical beneficial aspects of Raman spectroscopy is that a spectrum 
acquired using a laser source at 1064 nm will be the same as spectra acquired at other 
wavelengths (Note; the transmission function of different Raman systems vary, often 
resulting in Raman spectra having band intensity differences for identical samples 
acquired with different systems; in addition Raman differences may be due to the 
scattering dependence). While the Raman fingerprint band locations are nearly 
wavelength independent of the frequency of the excitation laser, the intensity of a normal 
Raman signal is inversely proportional fourth power of the frequency of the excitation 
source. This frequency dependence is evident in the following equation; 
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(3.1) 


/ = K/aV 


= K'la^ - 
A 


where I is the scattered signal intensity, K and K’ are composite constants, / is the laser 
power, a is the polarizability of the electrons in the analyte, and v and A are the 
frequency and wavelength of the incident light. For normal (non-resonance) Raman, the 
polarizability of the molecule is constant, but this term becomes very important for the 
resonance Raman processes. 

The first standoff Raman measurement of high explosives was obtained using a 
532 nm laser source [19]. It follows from taking the ratio of Equation 3.1 at two different 
wavelengths that shifting to a 229 nm laser source should achieve an immediate 29-fold 
improvement in the signal intensity. Raman spectroscopic measurements, acquired with 
excitation wavelengths in the near-infrared (NIR), are commonly performed to reduce 
contributions from sample fluorescence, which is more likely to occur when using 
excitation wavelengths in the visible portion of the electromagnetic spectrum. But a 
Raman measurement made at 229 nm should theoretically have a 138-fold advantage to 
that made using a 785 nm source, assuming all other factors are the same. However, 
Raman spectroscopy is not one-dimensional; other notable factors to consider include 
sample fluorescence, absorption, potential resonance enhancements, matrix interferences, 
and particularly laser-induced sample degradation. 

2. Resonance Raman 

In addition to frequency enhancements gained from using a UV laser to make 
Raman spectroscopy measurements, additional gains are achieved when the energy of the 
laser is close to the energy of a strong electronic transition of the analyte molecule. This 
enhancement is known as resonance Raman scattering and can provide enhancements of 
three to four orders of magnitude [18]. 
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Figure 6. Illustration of the prineipal difference between normal Raman and resonance 

Raman. The virtual state in normal Raman spectroscopy is not necessarily a true 
quantum state of the molecule, like those participating in resonance Raman. 

In theory, resonance Raman enhancements are most pronounced when the 
excitation wavelength is chosen to coincide with a wavelength () where the molar 
absorptivity of the analyte is maximum (); this is an electronic transition within the 
molecule. Resonance enhancements associated with a given electronic transition 
increases proportional to the square of the molar absorptivity [20]. 

Analysis of the Kramer Heisenberg Dirac (KHD) expression provides 
mathematical insight into the means of enhancing the resonance Raman effect [18]. The 
KHD expression is 

N ,^r (F|rJ/X/|rjG) </k,|G)(F|f. |/) ^ 

co„ + co,-ir, )’ 

where a is the molecular polarizability, p is the direction of the incident polarization, a 
is direction of the scattered polarization, G is the ground vibronic state, I is an excited 
vibronic state, F is final vibronic state of Raman scattering, is the dipole operator, 

is the frequency of the laser/incident radiation, ®g^and co^p are difference angular 
frequencies between the ground state and an excited vibronic state and the excited state 
and final state respectively, and iT j is small factor relating to the lifetime of the virtual 
excited state called the phenomenological damping constant [18]. Use of Dirac’s “bra- 
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ket” notation makes this summation over all vibronie states look deeeptively simple 
eonsidering that eaeh numerator eonsists of two triple integrals multiplied by one another. 
Fortunately this equation ean be easily disseeted [18]. The summation oeeurs over all 
vibronie states, although typieally only the first few exeited states are primarily involved. 
The numerators of the KHD expression mathematieally express the nature of the Raman 
seattering depieted in Figures 3 and 6. Sinee there is no reason why this proeess should 
have to start in the ground state, there are two terms required to aeeount for both Stokes 
and anti-Stokes seattering. The simplieity in analyzing this equation eomes from analysis 
of the denominator, eommonly ealled the resonanee denominator. The energy of the /F^ 
term, whieh is a term added to aeeount for the lifetime of the exeited state, is small 
eompared to co^j and . If the laser wavelength is ehosen to eorrespond elosely to the 
energy required to exeite the moleeule into an eleetronie transition, the denominator for 
Stokes seattering beeomes very small (approximately /F^, whieh itself is small), and the 
resulting value of the evaluated equation grows signifieantly aeeounting for the resonanee 
effeets. Not all Raman bands will be enhaneed over normal Raman proeesses [23]. The 
KHD equation is useful to gain physieal insight eoneerning the Raman resonanee effeet, 
but sinee the observed perturbations are large, quantum perturbation theory is not valid. 
For a detailed theory of Raman resonanee using a vibronie eoupling model see Ref [22]. 

Sinee most samples absorb in the UV wavelengths, self-absorption of the Raman 
seattering signal faetor eomplieates resonanee Raman measurements over normal Raman 
teehniques [21]. Additionally, this absorption ean eause photodegradation of the analyte. 
These effeets ean mitigate some of the potential benefits of resonanee Raman 
speetroseopy. 

C. NOISE 

1. Shot Noise 

Speetrometers that use a eooled CCD or intensified CCD (ICCD) are normally 
sample shot noise limited; the dark noise is normally very small in these deteetors [23]. 
Shot noise is a fundamental phenomenon that arises from the random nature of emission 
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of light in excitation processes. Shot noise will always be present in Raman 
measurements, but in standoff Raman measurements the background noise will typically 
dominate over the fundamental shot noise limits of the measurement. 

2. Background Effects 

Background noise is any photon detected by the system that does not emanate 
from the analyte. Background effects are typically signihcant noise contributors in 
Raman measurements. When concomitant species in a sample matrix have Raman shifts 
near the analyte peaks, peaks can get convolved and obscured. The simplest example of 
this is the atmospheric O 2 Raman band that can either convolve or strongly dominate over 
the band that arises from an asymmetric stretch R-NO 2 in solid explosive samples. 

When the excitation laser induces fluorescence of the analyte or in other species 
in the sample matrix, the very broad fluorescence peak can significantly raise the noise 
floor and in some cases completely obscure analyte peaks. Unlike the specific 
‘fingerprint’ Raman selection rules, fluorescence selection rules are less specific and 
typically do not provide as useful chemical information for qualitative identification of a 
molecule, particularly in mixtures of samples. If fluorescence effects are prevalent there 
are a number of techniques that can be employed including changing the excitation 
wavelength. Selecting a near infrared (NIR) laser or using a laser below 260 nm can 
significantly reduce or eliminate background fluorescence [21]. 

3. Signal-to-Noise Ratio 

The Raman signal intensity divided by the root-mean-square noise signal is the 
typical definition of the signal-to-noise ratio (SNR or S/N) [5]. This ratio is a key 
indicator of the quality of the Raman measurement, and typically the signal peak 
becomes visually imperceptible over noise signals when the SNR falls below a value of 
two or three [24]. In Raman spectroscopy applications where background signals are 
relatively consistent, the S/N ratio can typically be increased by increasing source 
irradiance, exciting the Raman process using lower frequencies of laser light, co-addition 
of a number of collected spectra, and using longer interrogation times [5, 24]. 


23 



THIS PAGE INTENTIONALLY LEET BLANK 


24 



IV. PREVIOUS RAMAN RESEARCH 


A. RAMAN DETECTION OF EXPLOSIVES 

1. Raman as a Forensic Tool 

Raman spectroscopy began to flourish in the 1970s with the increased availability 
of lasers. During this period, IR and Raman spectra of several explosives were published 
and the fundamental molecular vibrations were assigned to observed frequency bands 
[25, 26]. Later, in the 1990s, Raman spectroscopy found a new application in forensic 
science when it was applied to the detection of explosives [27, 28, 29]. While validated 
as a new analytical method to detect explosives, due to the unavailability of Raman 
spectrometers in most analytical laboratories in the 1990s, it became more widely used to 
confirm the presence of explosives after other analytical tests had been used to detect 
their presence. Other forensic applications grew out of these efforts including in-situ 
analysis with portable instruments and even the analysis of explosive residues on 
fingerprints [18]. 

2. Resonance Raman Spectroscopy of Explosives 

In the late-70s and early 80s, the application of resonance Raman in physical 
chemistry and analytical chemistry fields began to emerge [20]. Sands et al. first 
observed UV resonance Raman spectra of impure explosives and narcotics using the 244 
nm line from a frequency-doubled Argon ion (Ar^) laser-based Raman microscopy 
system [30]. Resonance enhancements greater than normal frequency (v"^) related 
enhancements were observed but not quantified. In some cases, simplifications of the 
Raman spectra were also observed as weaker bands were obscured. 

Nagli et al. determined the absolute Raman cross sections of TATP, PETN, RDX 

and TNT explosives from 620 to 248 nm at a constant flux of 2.5 • lO^"* quanta s ' cm^^ 

using KNO 3 (potassium nitrate) powder as an internal standard [31]. They reported 

normalized UV Raman signals were 100-200 times stronger at 248 nm when compared to 
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measurements made at 532 nm. Their study also ineluded reported strong fluoreseenee 
of many baekground materials at 532 and 785 nm exeitation with a signifieant reduetion 
in baekground fluoreseenee at 248 nm. 

Comaneseu et al. reeently published a eomprehensive study on the resonanee 
Raman speetra of TNT, RDX, HMX, and PETN using 40 UV wavelengths from 210 to 
280 nm using a 90° eolleetion geometry [32], This study ineludes tables that detail 
exeitation wavelengths that produeed the strongest Raman bands for HE samples in 
aeetonitrile solutions. The strongest Raman bands tend to be the symmetrie NO 2 streteh 
and the symmetrie phenyl streteh, whieh is eonsistent with Raman seleetion rules. 

B. STANDOFF RAMAN 

In 1992, Angel et al. demonstrated standoff Raman deteetion of ferroeyanide at 
distanees of 20 m using 488 nm and 809 nm eontinuous wave (CW) lasers, a 40 mm 
eolleetion optie, speetrograph, and a liquid-nitrogen eooled CCD [33]. 

In 2000, Ray et al. at Brookhaven National Eaboratory developed a UV Raman 
lidar system for the standoff deteetion of eontaminants at distanees of 30 m or less using 
a pulsed Nd;YAG laser operating at 266 nm with an average intensity of 200 mW/em , 
six ineh Newtonian teleseope, speetrograph, and CCD eamera [34]. In the artiele, Ray 
and eoworkers used eyelohexane, aeetonitrile, and Teflon® for his proof of eoneept. 

In 2005, Carter at Eawrenee Eivermore National Eaboratory, in eooperation with 
Angel and other researehers at the University of South Carolina, performed the first 
Raman standoff deteetion measurements of high explosives in a siliea matrix [19]. Using 
a pulsed 10 Hz Nd;YAG laser at 532 nm, an 8 in Sehmidt-Cassegrain teleseope, 
speetrograph and a gated ICCD eamera, various explosive samples were deteeted at 
standoff distanees up to 50 meters. This work also ineluded a power density study whieh 
demonstrated linear Raman intensities up to ~3 x 10^ W em'^. Additionally, in some HE 
samples photodeeomposition or other signal degrading effeets were observed. 

Sharma and eollaborators at the University of Hawaii published their standoff 
Raman results just a few months following Carter’s work. Sharma measured milligram 
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(mg) quantities of explosives samples at 10 m using a pulsed Nd:YAG laser (532 nm, 10 
Hz) with a Cassegrain-Maksukov teleseope (f/15) and ICCD [35]. 

Gaft and Nagli reeently reported sub-mierogram (sub- /ug) deteetion limits of 
explosives using a gated standoff Raman system at distanees up to 30 m [36]. 
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V. METHOD 


A. RAMAN CROSS SECTIONS 


The intensity of a Raman line ean be expressed in two different ways: (1) from a 
purely theoretieal or (2) an experimental viewpoint. In the former ease, the total Raman 
seattering eross seetion, a, for the vibrational transition from n<—m for an isolated 
moleeule deseribed in Ref [37] as follows: 

= (5.1) 


where Ima is the Raman intensity over 4n steradians of the vibrational transition n<—m 
divided by the ineident exeitation intensity, lo. The above expression is a funetion of the 
ineident and Raman mode frequeneies, the degeneraey and thermal oeeupaney of state m, 
and the Raman polarizability tensor [37]. In the latter ease, the experimentally observed 
Raman line intensity is given as the following expression adapted from Ref. [37]. 

L„=h-N-D-K-CT„ (5.2) 

where Imn is related to the moleeular density {N) in moleeules em' , ineident exeitation 
intensity (/o), sample pathlength {D - based on ineident beam and Raman seattered 
attenuation), several variables for eolleetion and deteetion of the seattered light {K), and 
the Raman seattering eross seetion {omn)- 


An absolute differential Raman eross seetion (cr^), a differential seattering 
probability per unit of solid angle, ean readily be ealeulated using internal standards and 
by eomparing the relative Raman eross seetion to a known absolute differential Raman 
eross seetion. This is shown in the following equation, whieh was adapted using 
aeetonitrile as the internal standard [37, 38]: 


{^rX ={(^rX,[J^ 


Vo(vo-919)^ 


I c. r- 


(5.3) 


29 



where (cr^)^ and are the absolute differential Raman eross seetions of the 

analyte and acetonitrile reference sample respectively, 1^^ and are the integrated 
scattering intensities of the analyte and acetonitrile bands respectively, Vq is the 
frequency of the incident excitation laser, 919 is the observed Raman shift of the 
reference acetonitrile peak in cm'\ is the frequency of the observed analyte band 
Raman shift in cm‘\ and Z is a dimensionless term that factors in the change in local 
field effects differences between gas samples and solutions [47, 48, 49]. L can be 
described using the index of refraction of the solution using the following expression [37, 
38,45,48,49]; 


L = 





(n^+lY 


(5.4) 


where the index of refraction of the solution (n ) is assumed to approximately equal the 
indices of refraction evaluated at the analyte band frequency (\-y ) and the frequency 

of the incident laser (n ). 

B. RAMAN ENHANCEMENTS IN SOLUTION 

Acetonitrile was chosen as an internal standard for this enhancement study for 
two reasons: it is an effective solvent for high explosives and it is a useful internal 
standard to quantify resonance enhancements. The 919 cm'^ C-C stretch in acetonitrile 
was shown by Dudik et al. [37] to be solely dependent on the well-known Raman 
frequency dependence since it does not display any resonance or pre-resonance 
enhancements. For this reason, Raman enhancements can be readily calculated by 
making Raman cross section calculations using the internal standard. 

To develop an expression for resonance enhancements. Equation 5.3 is evaluated 
at pre-resonance/resonance Raman conditions (UV excitation denoted by subscript UV) 
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and divided by the same equation evaluated under normal (i.e., non-resonanee) Raman 
eonditions (NIR or visible exeitation with the subseript NIR). This yields: 
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Sinee approximately varies only with thev"^ frequeney enhaneements, 


[(^i^)919_ 


^ur 

[(^i^)919_ 

'^NIR 

_ ^NIR _ 


(5.5) 


(5.6) 


Equation 5.5 ean be rearranged to give an expression for the resonanee Raman 
enhaneement of UV measurement relative to the normal Raman measurement that is 
independent of the v‘^ Raman frequeney enhaneements: 
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(5.7) 


This step assumes that the index of refraetion of the solution used in the NIR 
measurements is approximately equal to the index of refraetion used in the UV 
measurements to eaneel the loeal field terms (T). The index of refraetion ean be 
expressed as a funetion of eoneentration and wavelength [39]. The index of refraetion is 
only weakly wavelength dependent and the eoneentration of all samples is approximately 
19 mol L"' (or M) aeetonitrile, so any error introdueed from the loeal field term will be 
small. Additionally, sinee the eoneentration of aeetonitrile is nearly the same for UV and 
NIR measurements, the terms eaneelled. 
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The Beer-Lambert law describes the transmission of light through a sample: 


T = 



(5.8) 


where T is the transmission, 1^ and / are the light intensities before and after entering the 
sample, s is the molar absorptivity (L mol'^cm’'), b is the path length (cm), and c is the 
concentration (mol L'^) [5]. The transmission of the sample may also vary across the 
Raman spectrum, and this effect is more pronounced in UVRR analysis. To derive 
Equation 5.3, Dudik et al. assumed that the transmission of the solution at the analyte and 
reference bands was approximately the same; however, it is necessary to reintroduce the 
transmission term since most explosive samples strongly absorb UV radiation. This leads 
to the final version of the Raman enhancement formula for solutions (E^^^^): 





“(Vo-919)^' 
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After Raman and EiV-Visible (ETV-Vis) spectra are collected and analyzed. 
Equation 5.9 is used to calculate the excitation wavelength dependent resonance Raman 
enhancements for individual Raman bands. 


C. RAMAN ENHANCEMENTS IN SOLIDS 

Since many high explosive solids highly absorb UV light, their molar absorptivity 
can be over lO"^ E mof'cm'^ When the concentration, molar absorptivity, and 
concentration of a sample is fixed. Equation 5.8 can be used to calculate the maximum 
penetration depth of UV laser light into a sample. This depth is typically of the order of a 
tenth of a micron to microns for explosive samples in UVRR measurements. If the 
photon is Raman scattered before being absorbed by the sample, the scattered Raman 
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signal can be self-absorbed before getting out of the sample (see Figure 7). Using 
Equation 5.8 and assuming no atmospheric attenuation, the strength of the laser at the 
point of Raman scattering is 

/ = (5.10) 

where Sq is the molar absorptivity of the sample at the wavelength of the incident laser, 
and b is the path length through the solid to the seatterer. Applying the Beer-Lambert 
law to the Raman signal out of the solid sample () at a particular path length, we now 
have 






= /o -e 


-2.3036c(£-o ) 


(5.11) 


where is the molar absorptivity of the sample at the wavelength of Raman signal and 

the path length from the seatterer back out of the sample is for a 180° collection 
geometry. 



Figure 7. Difficulties arise in UVRR due to the attenuation of the laser into the sample 
and the self-absorption of the Raman signal out of the sample. 


The total magnitude of Raman signal over all path lengths for a given sample 


depth (D) is also proportional to the differential eross seetion of the molecule (cr^ ) in 
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2 1 3 

cm sr' (steradian), molecular density (N) in molecules em' , and K deseribes the 
variables for oolleetion and deteetion of the seattered light [40]; 


I 




(5.12) 


The physical interpretation of the integral in Equation 5.12 is the penetration 
depth of the sample, where Imn is reduced to a value of e'^ (-0.368), whieh only equals the 
sample depth when molecular absorptivity is negligible [40]. Taking the ratio of UV and 
NIR measurements using Equations 5.9 and 5.12 is a step closer to deriving a formula for 
Raman enhaneements in solids (E^): 




)o 

C/E 

Ivffi 
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1 ^ ^-2.303Z»c(£-o )^^ 
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NIR 


Additionally, it follows from Equations 5.6 and 5.7 that; 




(5.13) 


(5.14) 


This leads to the final relationship for eomparing observed Raman enhancements 
in solution to solid samples; 


^SOLN 



-2.303Z)c(fo+£'^)^^ 


-iUV 


-2.303dc(£'Q ) 


db 


NIR 


(5.15) 


where E^and are the local field effeet terms evaluated for solid and solution 

samples respeetively. In Equation 5.14, eq and 8 r represent the eombined absorption and 
scattering eoeffieient for the solid (e.g., powder or polyerystalline) sample. In this final 
step a correetion term for E, whieh is expeeted to be stronger in solids than in solutions, 
was applied. E^ and E^^^^ can be ealeulated using the wavelength dependent indices of 
refraetion, some of whieh are available in the literature, e.g., see Ref [40]. Using 
Equation 5.14, the resonance Raman enhancements from solution measurements can now 

be analyzed to estimate the magnitude of total enhaneements (i.e., resonanee Raman 
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enhancement multiplied by the frequency dependent enhancement and local field 
enhancement) of UVRR measurements in solid samples. 

Finally, if these measurements are to be applied to standoff distances, the 
atmospheric attenuation of the laser and Raman signal in air over a standoff distance 
(Standoff) i^ust be considered: 

-2rf^£,C; 

"^Standoff ^ ^ mn ^ ^ ‘ ? (5.16) 

where d is the standoff distance and s. and Ci are the molar absorptivity and 
concentrations of individual components of the atmosphere respectively. For UV 
wavelengths, ozone (O 3 ) has been reported as the principle absorbing species [14]. 


35 



THIS PAGE INTENTIONALLY LEET BLANK 


36 



VI. EXPERIMENTAL TECHNIQUES 


A. OBJECTIVES 

The goal of this research was to quantitatively determine resonance Raman 
enhancements in the absence of the frequency dependent enhancements for select high 
explosive (HE) materials by contrasting Raman measurements at 785, 244, and 229 nm 
excitation wavelengths using an internal standard. The internal standard (acetonitrile) 
exhibits no resonance or pre-resonance Raman effects and did not overly convolve the 
analyte Raman peaks. The experimental conditions between ultraviolet resonance Raman 
(UVRR) and near-infrared (NIR) Raman systems were matched to the greatest extent 
possible given experimental constraints, and the use of an internal standard allows 
quantitative comparison between measurements. Using this technique, it is possible to 
determine the benefits of UV excited Raman compared to normal Raman techniques for 
the standoff detection of high explosives. 

B. EXPERIMENTAL SUMMARY 

UV-visible (UV-Vis) absorption measurements were performed to select solution 
concentrations for UVRR analysis and to collect required transmission data for Raman 
enhancement calculations and corrections. Sample solutions at a single fixed dilute 
solution (varying from 2.2 to 9.27 mM for different analytes) were placed in a 1.00 mm 
path length (sealable) quartz cell for UVRR analysis. UVRR measurements of sample 
solutions were made for three to five sample replicates at five different laser powers 
settings for each replicate ranging from 20 juW to 2.2 mW. The laser was rastered 
across the cell during the 180-second analysis time to prevent or minimize localized 
heating and photodegradation effects that could negatively influence the measurement. 
Three replicates of one fixed sample concentration between 100 to 150 mM were used for 
NIR Raman measurements at 50 mW. The higher concentration was necessary since the 
analyte Raman band intensities were weaker for NIR excitation; NIR excitation provides 
no resonance enhancements effects and Raman scattering power goes as v^. The laser 
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power for NIR measurements did not produee any observable degradation of any 
samples; Raman intensity measurements over six minutes for fixed spot were found to be 
highly reprodueible. 

Approximately 10 mg quantities were used for Raman analysis of solid samples. 
Samples were plaeed on glass slides for UVRR and NIR Raman analysis for 180 seeonds. 
The UV laser was rastered aeross the solid samples for UVRR analysis to prevent 
photodegradation effeets and heating. Raman speetra were analyzed to determine 
integrated peak intensities (area under the Raman bands) to perform enhaneement 
ealeulations. 

The use of the internal standard greatly simplifies the eomparison of 
measurements made on separate Raman systems. Taking the ratio of the analyte and 
referenee bands in Equation 5.8 negates the differing power settings between the systems, 
differing sampling volumes, differing deteetor pixel sizes, differing system f-numbers, 
sample absorbanee of the exeitation laser in the sample, and other differenees that eannot 
easily eontrolled experimentally. The speetral bandpass of the two Raman systems and 
sample analysis times were matehed for added measure. 

C. EQUIPMENT 

Two separate Raman systems were used at Lawrenee Livermore National 
Laboratory to eolleet the data for this thesis. Measurements were eolleeted using 
ultraviolet and near infrared exeitation for resonanee and normal Raman speetroseopy. 
All laser power measurements were made using an Ophir NOVA II power meter 
(Bioneties, St. Laurent, QC, Canada). 

Sample absorbanees were measured using a Cary 300 Bio UV-Vis 
speetrophotometer (Varian, Ine., Palo Alto, CA, USA) and baseline eorreeted using 
aeetonitrile. 
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1 . 


NIR (785 nm) Raman System 


This Raman system consisted of an Invictus 785-nm NIR Laser (Kaiser Optical 
Systems, Ann Arbor, MI, USA), Mark II Filtered Probe Head utilizing holographic optics 
(Kaiser Optical Systems, Ann Arbor, MI, USA), f/4 spectrograph (Chromex, Model 
250IS/RF, Albuquerque, NM, USA), liquid-nitrogen-cooled CCD system (1340x400 
pixels, Princeton Instruments, Trenton NJ, USA). The Raman signal intensity captured in 
the vertical channels of the ICCD was binned to a single channel. The spectrometer and 
optical fiber probe had focal ratios of f/4 and f 12, respectively, resulting in some Fresnel 
losses. The spectral bandpass of the spectrometer was 12 cm'' for a 149 /rni slit width 
and the 600 grooves / mm grating. 



Figure 8. Picture of the NIR Raman System showing the placement of key components. 

This system was calibrated using a 50% v/v (volume / volume) solution of toluene 
/ acetonitrile (see Appendix A). Peaks locations on the camera’s pixels were matched to 
known Raman shift locations in wavenumbers for the reference standard using a 
quadratic peak-fitting model using WinSpec 32 software. 
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2 . 


UV Raman Microscopy System 


This device consists of a water-cooled Argon ion laser (Lexel, Model 85-SHG, 
Fremont, CA, USA), microscope (Leica DMLM, Leica Microsystems, Wetzlar, 
Germany), ICCD (2048 x 512 pixels, Princeton Instruments, Trenton, NJ, USA), and 
f/5.7 spectrograph (Model PI SpectraPro 25001, Princeton Instruments, Trenton NJ, 
USA). The spectral bandpass of the spectrometer was 12 cm'^ for a 100 /um spectrometer 
slit width and the 2400 grating / mm grating. The microscope translation stage was 
equipped with an electronic joystick controller for manual x- and y-axis rastering of the 
laser across the sample. 



Figure 9. Photo of the UV Raman Microscopy System highlighting the system’s key 

components. 

The UV laser system was initially setup for 244 nm (frequency doubled 488.0 nm 
line) excitation with a maximum power at the sample of 2.2 mW. The laser was later 
modified to emit at 229 nm by frequency doubling the 457.9 nm line. The addition of 
new optics suitable for UV wavelengths, including the addition of two Pellin-Broca 
prisms to remove Ar ion plasma lines, reduced usable power to 1.0 mW at the sample for 
229 nm measurements. 
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The manufacturer's wavelength calibration was used for all UVRR system 
measurements at 244nm excitation. For 229 nm measurements it was necessary to 
calibrate using cyclohexane as a Raman shift standard (see Appendix A). 

D. SAMPLES AND SAMPLE PREPARATION 

Toluene from EMD Chemicals Inc. (Darmstadt, Germany) and spectroscopic 
grade acetonitrile from Riedel-de Haen Laborchemikalien GmbH & Co. (Hannover, 
Germany) was used as received. Acetonitrile was used for all Raman solution 
preparations. Cyclohexane for select UV-Vis spectrometer measurements was obtained 
from Sigma Aldrich (St. Louis, MO, USA). 

2,4-DNT (2,4-dinitrotoluene); 2,3-DNT (2,3-dmitrotoluene); 3,4-DNT (3,4- 
dinitrotoluene); 2,6-DNT (2,6-dmitrotoluene); PNT (4-nitrotoluene), 2,4,6-TNT (2,4,6- 
trinitrotoluene), HMX (l,3,5,7-tetranitro-l,3,5,7-tetrazocane), RDX (1,3,5- 
trinitroperhydro-l,3,5-triazme). Composition B (63% RDX, 36% TNT and 1% paraffin 
wax), PETN (pentaerythritol tetranitrate), ANEO (ammonium nitrate / fuel oil), and 
TATB (l,3,5-triammo-2,4,6-trmitrobenzene) were high purity samples provided by 
LLNL’s High Explosives Applications Eacility (HEAE) and Eorensic Science Center. 

All reagents, with the exception PETN, were used without further purification. 
PETN solutions were centrifuged at 6000 rpm for 12 minutes to remove scattering 
particulates found in 150 and 9 mM solutions. PETN is typically coated with 
polyethylene to reduce the sensitivity of the explosive, and UV Raman analysis of the 
PETN impurity insoluble in acetonitrile was consistent with polyethylene and other 
similar molecules. 

Solid samples (approximately 10 mg quantities) were transferred from storage 
vials to glass slides (Gold Seal, Thermo Eisher Scientific Inc., Waltham, MA, USA) for 
spectral analysis. 

The same UV/Vis cuvette (1.00 mm path length, QS grade, Hellma, Eorest Hills, 
NY, USA) was used for all solution spectral measurements. The cuvette was triple rinsed 
with acetonitrile, air dried, and rinsed with the new sample prior to introducing a new 
sample for measurement. 
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E. DATA HANDLING 

WinSpec 32 spectroscopy software (Princeton Instruments, Trenton, NJ, USA) 
was used for data acquisition and system calibration. Data fdes were converted to ASCII 
files and plotted using IGOR Pro 6 (WaveMetrics, Portland, OR, USA). Data from UV 
Raman measurements were recorded in wavelength (nm) and converted to Raman shift in 
wavenumbers (cm'^) using a simple MATLAB (The Math Works, Natick, MA, USA) 
program. 

The goal of making Raman band area calculations was to use an accurate and 
simple program that could be used by most analytical laboratories without requiring 
special software. MATLAB programs (see Appendix B for source code) were used to 
baseline the spectra between two user-selected points, to select the start/end of a Raman 
peak, and to calculate the area under spectral peaks for subsequent Raman enhancement 
calculations. Two versions of the program were utilized; one to be used with 
deconvolved peaks and another that can be used with peaks that were partially convolved. 
The latter program calculated half of the peak area of a visible peak and invoked 
symmetry to calculate the full peak area of the band; a similar method has been reported 
by Mallard et al. [46]. The use of PeakFit (Aspire Software International, Ashburn, VA, 
USA) nonlinear curve fitting software was also summarily examined for baselining and 
deconvolution of the data; however, the Gaussian, Lorentzian and other user-selected 
data models did not accurately fit the form of the peaks obtained in this research. The 
Raman peaks should be mixed Gaussian and Lorentzian functions due to pressure 
broadening, Doppler sifting, and other instrument response features [41]. PeakFit may 
still be an effective tool for these calculations, but this research emphasized the 
developed of MATLAB tools and so the use of Peakfit was abandoned. The accuracy of 
the MATLAB model was validated using the cumulative area tool in PeakFit, which does 
not require fitting the peaks to any mathematical models. 
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VII. RESULTS 


A. SUMMARY OF RESULTS 

1. Molar Absorptivity and Solvent Effects 

Resonance Raman enhancements of PETN; PNT; 2,4-DNT; and 2,4,6-TNT were 
measured using acetonitrile as both the solvent and internal standard. In theory, 
resonance Raman enhancements are most pronounced when the excitation wavelength is 
chosen to coincide with a wavelength ( ) where the molar absorptivity of the analyte 

is maximum this is an electronic transition within the molecule. Resonance 

enhancement associated with a given electronic transition increases with the square of the 
molar absorptivity [20]. The molar absorptivity of high explosive samples and select 
explosive precursors dissolved in acetonitrile were calculated from absorption 
measurements and are presented in Figure 10. 



Figure 10. Calculated molar absorptivities (s) from absorption measurements of PFTN, 
PNT, 2,4-DNT, and 2,4,6-TNT in acetonitrile. It is interesting to note that the 
molar absorptivity 2,4,6-TNT > 2,4-DNT > PNT at which corresponds with 
the increasing number of NO 2 substitutions on the aromatic ring. 
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The observed results for the nitroaromatie speeies were eompared to 


literature values and are presented in Table 2. 


Analyte 

Literature 

Observed 

4ax (nm) (M 'em ') 

4ax (nm) (M 'em ') 

PNT 

2,4-DNT 

2,4,6-TNT 

264 10100 

233 15800 

224.5 23000 

276 9470 

245 13900 

217.5 20900 


Table 2. Caleulated and referenee maximum molar absorptivities for the primary n ^ n* 
transitions observed at wavelengths >220 nm for PNT; 2,4-DNT; and 2,4,6-TNT 
in aeetonitrile. Referenee values for PNT and 2,4-DNT were obtained using n- 
heptane as a solvent and 2,4,6-TNT was observed in eyelohexane [42], 

The results obtained for and for the n ^ n transition (a stable 
bonding to an unstable anti-bonding orbital) for PNT and 2,4-DNT notieeably differ from 
the literature values due to solvent effeets. These analytes display a bathoehromie (red 
shifting) of the UV-Vis speetra in the polar aeetonitrile solvent. The observed 
bathoehromie shifts are eonsistent with the polar aeetonitrile solvent moleeules 
stabilizing (lowering) the exeited state relative to the ground state [43]. The 
bathoehromie shift is very notieeable in the following example experimental results of 
ealeulated molar absorptivities for 2,4-DNT in polar aeetonitrile and non-polar 
eyelohexane solutions (Figure 11). 


44 






Wavelength (nm) 

Figure 11. Calculated molar absorptivities (£■) of 2,4-DNT in acetonitrile and 
cyclohexane. The use of a polar solvent produces broadening and bathochromic 
shifting of the peak of the main;r ^ n transition observed at 234 nm in 
cyclohexane and 245 nm in acetonitrile. Also note that some hypochromicity 
(absorption decrease) is noted for the polar solvent. 


It is more difficult to determine if a frequency shift occurs with 2,4,6-TNT due to 
its broad and relatively flat peak maximum, but the data may suggest a hypsochromic 
shift (blue shift) of the n ^ n transitions, which can also be observed in polar solvents 
where the ground state is lowered relative to the excited state [43]. Due to the observed 
bathochromic and possible hypsochromic shifts of n ^ n electronic transitions of 
analytes in acetonitrile solutions, the wavelength for maximum Raman resonance 
enhancements may differ for the same HE materials in the solid phase. For this reason, 
previously observed resonance maxima by Conanescu et al. [32] may differ slightly 
(perhaps up to 10 nm) from solid samples. Additionally, the calculated values of Raman 
enhancements at particular wavelengths presented in this thesis may also differ slightly 
from actual values due to the same solvent effects shifting . 

2 . Solid Explosive Raman Spectra 

Raman spectra of select high explosive samples are shown in Figure 12 using a 
785 nm excitation laser. Additional spectra can be found in Appendix A. Due to specific 
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vibrations that are influenced by the local chemical environment, each Raman spectrum 
is generally unique (i.e., a “fingerprint”). Military explosives are most commonly used in 
formulations with mixtures of other explosives and compounds to improve mechanical 
and chemical properties [32], In Figure 12, the Raman spectral signature of Composition 
B, a mixture of TNT and RDX, is plotted along with pure samples of RDX and TNT. 
Although at first glance the Composition B Raman spectrum looks like that of TNT, a 
closer inspection reveals the presence of RDX peaks at ~890, -1270, and -1310 cm"'. 
Composition B’s Raman spectra also exhibits some broadband fluorescence which is 
likely due to the 1% paraffin wax it contains. 



Figure 12. Raman spectra of RDX, TNT and Composition B explosives acquired using a 
785 nm laser at 50 mW. TNT and RDX are the principal explosive ingredients of 
Composition B, and their composite signatures are visible in the Composition B 
spectrum. TNT and RDX spectra are offset for clarity. 

When Raman spectra are collected using an excitation wavelength at or near a 
strong electronic transition of the target molecule, the resulting spectrum is often 
different (even simplified) compared to a spectrum acquired with wavelengths far from 
strong transitions. The differences between these spectra are (in the absence of 
degradation) a result of resonance enhancement in which only the molecular vibrations 
associated with the electronic transition are amplified and appear in the acquired Raman 

spectrum. The other analyte peaks may be present but not visible if the enhancement is 
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strong. Since most molecules absorb in the ultraviolet spectral region, resonance Raman 
effects are often observed when using UV excitation wavelengths. For example, some 
Raman spectral simplification can be observed in the ultraviolet resonance Raman 
(UVRR) spectra of TNT and PETN compared to normal Raman as shown in Figures 13 
and 14. This simplification can make qualitative determination easier in complicated 
sample matrices, but it can also make it more difficult when trying to distinguish between 
molecules with very similar molecular structures. 



Raman Shift (cm 

Figure 13. Raman spectra of PETN using 785, 244, and 229 nm laser excitation 

wavelengths. 


47 













Raman Shift (cm ’) 


Figure 14. Normalized Raman speetra of 2,4,6-TNT using a 785, 244, and 229 nm laser 
excitation wavelengths. The asymmetric NO 2 peak is clearly visible in the NIR 
Raman spectrum, but the atmospheric oxygen band (~1585cm'^) in the 229 nm- 
excited spectrum fully obscures this analyte band. 

The presence of the atmospheric O 2 Raman band in the UVRR spectrum is 
problematic for observing the asymmetric NO 2 stretch in most of the samples used in this 
study. For real-world standoff measurements, the O 2 band will most likely be the 
dominant Raman band in this spectral region and completely mask the weaker 
asymmetric NO 2 Raman band of the nitroaromatic explosives. 

Figure 15 demonstrates that with sufficient detector resolution and S/N it is 
possible to unambiguously identify analytes with very similar molecular structures using 
UVRR. 2,4-dinitrotoluene and 2,6-dinitrotoluene only differ by the location of one NO 2 
group on the molecule’s aromatic ring; still these spectra can be differentiated with 
sufficient S/N. 
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Figure 15. Raman spectra of PNT; 2,3-DNT; 2,4-DNT; 2,6-DNT; and 3,4-DNT crystals 
collected using a 229 nm Raman microscope. Laser power was 0.27 mW with 
180-second integration time for all spectra. Spectra are offset by 2000 counts for 
clarity. Sharp atmospheric O 2 and N 2 Raman lines are indicated. 

Due to resonance effects, it is also possible to observe additional analyte 
vibrational bands in UVRR compared to normal Raman by amplification of bands that 
are otherwise very weak. This is a likely explanation for the appearance of additional 
peaks observed in ANFO’s UVRR spectra shown in Figure 16. It is also common to 
observe overtones of other bands in UVRR spectra. 
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Raman Shift (cm') 

Figure 16. Normalized Raman speetra of ANFO using 785, 244, and 229 nm laser 
excitation wavelengths. The Raman bands between 1200 to 1700 cm’' are more 
pronounced at lower excitation wavelengths. 

The acquisition of Raman spectra for RDX proved to be very challenging for UV 
excitation, even at low laser power settings while rastering the sample. The UVRR 
spectra of solid RDX, a molecule that is less stable than TNT, exhibits various broadband 
features that do not map well to any of the Raman spectral features in the NIR spectrum, 
as shown in Figure 17. The very poor results obtained for RDX in the UV may be a 
result of photodegradation of the molecule. Our UV results suggest that RDX would 
likely be difficult to detect in a sample matrix due to the rather ambiguous observable 
spectral features. However, RDX is most commonly used in formulations with other 
explosives (Comp B and C-4 are notable examples), so this drawback may not be a 
hindrance to tactical explosive detection efforts. 
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Raman Shift (cm’^) 

Figure 17. Raman spectra of RDX using 785 nm, 244 nm, and 229 nm laser sources. The 

peak near 1550 cm"' is atmospheric O 2 . 

3. Raman Spectra of HE Solutions 

Raman intensity measurements of the same sample while cycling through two 
laser power cycles from 0.025 to 2.2 mW were highly reproducible (<2% error). Raman 
measurements with the laser pointed (i.e., no rastering of the sample) at a fixed spot in 
PNT solution demonstrated an approximate 18% reduction in Raman intensity after live 
180-second scans. This reduction in signal is likely due to index of refraction changes 
from localized sample heating and/or photodegradation. Rastering the laser across the 
solution samples produced highly reproducible Raman spectra for all samples and was 
the technique utilized for all quantitative measurements. 

Since solute concentrations ranged from only 2.2 mM for UVRR analysis to 150 
mM for NIR Raman analysis, the solvent (acetonitrile) Raman bands expectedly 
dominated the observed Raman spectra for sample solutions. Raman signatures of neat 
acetonitrile and PETN solution in acetonitrile are presented in Figure 18. The 
asymmetric and symmetric NO 2 peaks of PETN are clearly visible in the PETN solution 
spectrum. 
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Figure 18. Raman spectra of 9.25 mM PETN solution in acetonitrile (bottom) and neat 
acetonitrile (top) at 244 nm excitation. Laser power was 100 pW and the sample 
analysis time was 180 s. The 1282 cm"' and 1650 cm'* bands of PETN are small 
but visible and marked by an asterisk. 


The UVRR spectra of TNT acquired at 244 and 229 nm excitation wavelengths 
are presented in Eigures 19 & 20, respectively. The symmetric NO 2 stretch (~1340cm'*) 
for the all of the aromatic compounds (i.e., PNT, DNT, TNT) studied were partially 
convolved with an acetonitrile peak. 



Raman Shift (cm’’) 

Eigure 19. Raman spectra of 2.2 mM 2,4,6-TNT in acetonitrile recorded at various laser 
power settings ranging from 25 pW (bottom) to 2.2 mW (top) using 244 nm laser 
excitation. Total accumulated CCD integration time for all spectra was 180 
seconds to allow quantitative analysis of the integrated Raman band intensities. 
Symmetric and asymmetric NO 2 stretches and phenyl ring stretches of TNT are 

annotated. 
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Figure 20. Raman spectra of 2.2 mM 2,4,6-TNT in acetonitrile recorded at various laser 
power settings ranging from 25 pW to 0.95 mW using 229 nm laser excitation. 

The insets are an expanded view of the NO 2 and phenyl ring Raman modes 
regions. Total accumulated CCD integration time for all spectra was 180 seconds. 

The integrated intensities for the TNT NO 2 and phenyl ring bands and the 919 
cm"' acetonitrile (internal standard) band were measured. The Raman intensities of the 
919 cm"' reference peaks were all linear with increasing power as predicted by Equation 
3 . 1 . The regression coefficient (R ) was greater than 98%. The regression coefficients 
for PNT and 2,4-DNT analyte peaks were also highly linear. However, the regression 
coefficient of Raman band intensity versus power for the TNT NO 2 symmetric band is 
nonlinear with increasing laser power (see Figure 21). 
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Power (mW) Power (mW) 

Figure 21. (a) Integrated intensity of the 919 em'^ aeetonitrile referenee band vs. laser 

power at 244 nm exeitation. The result is highly linear as demonstrated by the 
regression eoeffleient indieated on the graph, (b) Integrated intensity of the 2,4,6- 
TNT symmetrie NO 2 streteh vs. laser power for 244 nm exeitation. This result is 
visibly nonlinear and ean be fitted with an exponential funetion 
(eonstant - exp(-/)). 


The nonlinear results of integrated line intensities versus power ean likely be 
attributed to photodegradation of the sample, but eould theoretieally be due to other 
nonlinear effeets that will be discussed later. 

B. CALCULATED RAMAN ENHANCEMENT 

The integrated Raman band intensity data combined with UV-Vis transmission 
data were used to calculate (refer to Equation 5.8) resonance Raman enhancements for 
various analyte bands in the Raman spectra acquired at 244 and 229 nm excitation (see 
Figures 22 to 25). These reported Raman enhancements are solely due to resonance 
effects. It is important to note that the additional multiplicative enhancement gained from 
the well-known frequency dependence (the integrated Raman scattering goes as v"^) of the 
Raman scattered signal have been factored out in Equation 5.8 by ratioing to the 
acetonitrile reference bands. The Raman enhancements for explosive analytes were 
generally found to decrease though a weak exponential function with increasing laser 
power. This is clearly demonstrated in Figure 22 with the phenyl ring stretch of TNT. 
This effect is most likely attributed to photodegradation, but could also arise from 
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resonance dampening or saturation. If photodegradation competes with the Raman 
scattering process, the Raman scattering cross section will be reduced and enhancements 
will be decreased. Saturation effects are normally associated with high-pulsed power 
lasers and depletion of the ground state. A shift in the center frequency of the resonance 
curve at higher powers due to other nonlinear effects (such as nonlinear damping) is also 
possible but not likely in this case. 

The uncertainty associated with the Raman resonance enhancement calculations 
using Equation 5.8 is estimated to be 16-18% based on propagation of experimental 
errors. The largest relative standard deviations are from Raman band area calculations 
and transmission measurements in the NIR region. 




Excitation Power (mW) 


Figure 22. Calculated Raman resonance enhancements for 2,4,6-TNT symmetric NO 2 
stretch (left) and phenyl ring stretch (right) at 229 nm and 244 nm laser excitation 

using Equation 5.8. 


For TNT, the most notable difference in resonance enhancements between the two 
UV wavelengths is with the Raman signal corresponding to the phenyl ring stretch. This 
stretch corresponds to an electronic transition centered at 260 nm in the UV-Vis spectra, 
so it is reasonable that enhancements at 244 nm are greater than those at 229 nm. The 
Resonance enhancements calculated for TNT’s symmetric NO 2 stretch are 
indistinguishable due to experimental error. 
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Figure 23. Calculated Raman resonance enhancements for PETN symmetric NO 2 stretch 
(left) and asymmetric NO 2 stretch (right) at 229 nm and 244 nm. 

For PETN, it can be predicted from the UV-Vis spectrum in Eigure 10 that there 
will be greater enhancements at 229 nm excitation compared to 244 nm for PETN, which 
was also observed experimentally for both the symmetric and asymmetric NO 2 stretch 
(see Eigure 23). 

The resonance enhancements for 2,4-DNT’s phenyl ring stretch at 229 nm and 
244 nm are generally indistinguishable owing to experimental error (see Eigure 24). 
Based on the molar absorptivity of 2,4-DNT, it is predicted that 244 nm excitation will 
produce larger Raman resonance enhancements in the phenyl ring band than at 229 nm. 
The Raman enhancements for the NO 2 stretch were greater at 229 nm than at 244 nm. 
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Figure 24. Calculated Raman resonance enhancements for 2,4-DNT at 229 nm (left) and 

244 nm (right). 


The molar absorptivities on PNT at 229 nm and 244 nm are close in value since 
there is a minimum in the UV-Vis spectrum at 235 nm; however, excitation at 244 nm is 
closer to the n ^ n transition associated with the phenyl ring stretch. As a result, the 
observed Raman enhancements for the PNT phenyl stretch is more enhanced at 244 nm 
than 229 nm excitation. 
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Figure 25. Calculated Raman resonance enhancements of symmetric NO 2 stretch (left) 
and phenyl ring stretch (right) for PNT at 229 nm and 244 nm laser excitation. 
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Using the method presented in this thesis, it is also possible to qualitatively 
demonstrate the Raman resonanee enhaneement eurves for 2,4,6-TNT as a funetion of 
wavelength by analyzing the data presented by Comaneseu et al. [32], The ratio of TNT 
Raman signal intensity to the aeetonitrile (919 em'*) ean be ealeulated direetly from their 
presented data, but there is not enough information to aeeurately ealeulate quantitative 
resonanee eurves. Figure 26 demonstrates the relatively large bandwidth of the Raman 
resonanee effeet for TNT. 




Figure 26. Caleulated Raman resonance enhancement curves for 2,4,6-TNT’s symmetric 
NO 2 band (left) and phenyl ring stretch (right) using acetonitrile as an internal 
standard. Data from; [32]. Relative analyte / reference peak heights were fitted 
using a Gaussian curve. The wavelength of maximum Raman resonance effects 
for TNT’s symmetric NO 2 stretch and phenyl ring stretch were reported to be 263 

nm and 260 nm respectively [32]. 


C. APPLICATION TO SOLIDS 

There were several reasons why it was desirable to demonstrate a method to 
estimate Raman enhancements in solution and apply those results to solid phase samples. 
High explosives in solution are much easier to deal with than solid samples for safety 
reasons and many labs, including LLNL, have limitations on the quantities of explosives 
that can be used in certain facilities. Additionally, resonance enhancement studies 
require an internal standard reference, which is much easier to introduce in solution- 
based studies; internal standards can be dissolved in the solution or the solvent can serve 
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as the internal standard. For solids, the introduction of an internal standard is more 
difficult since heterogeneities in terms of particle sizes and distribution are present. This 
results in irreproducible measurements and large error especially when one considers that 
the typical laser beam spot size is on the order of tens of microns in diameter. 

The effects of laser attenuation and self-absorption of the Raman signal were 
derived in the Methods section. The laser penetration depth and return Raman signal are 
functions of the optical (absorption and scattering coefficients) and physical (particle size 
and sample thickness) properties of the sample. The wavelength dependent laser 
penetration depths, evaluated using the integral inside Equation 5.14, are shown in Figure 
27 for 2,4,6-TNT and PETN. Since the wavelength of the Raman signal of the analyte 
varies with the excitation wavelength, the value of the molar absorptivity, SJ^, at each 
desired excitation wavelength must first be determined. Note: The biggest u nk nown 
material property, and the one most difficult to characterize, is the light scattering 
coefficient, which was not included in the model of Equation 5.14. It is recognized that 
future model improvements should include the scattering coefficient. 




Figure 27. Calculated penetration depth of the UV Raman signal for the phenyl ring 
stretch of 2,4,6-TNT (left) and symmetric NO 2 stretch of PETN (right). 

From Figure 27, it is clear that the deep-UV excited Raman scattering from 
explosives will come from very near the sample surface. In general, UV excitation 
provides a higher Raman cross section for a sample, because of the dependence of the 
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Raman intensity on the inverse 4th power of the wavelength. However, for highly 
absorbing samples, or sample matriees, the exponential deerease of the penetration depth 
of the laser at UV wavelengths leads to an overall deerease of the measured signal. It 
may seem eounterintuitive but in eases where resonanee Raman enhaneements are also 
realized, the optimal wavelength seleetion for maximizing the resonanee enhaneement 
may not neeessarily be the optimal wavelength for aehieving the highest Raman signal 
and sensitivity from the sample. As an example, NIR exeitation has an obvious 
advantage over UVRR for non-absorbing bulk samples (i.e., thiek samples) sinee this 
wavelength will penetrate deeper into the sample and provide more seattering that ean be 
deteeted. 


The resonanee Raman enhaneement ealeulation, using Equation 5.14, for the TNT 
phenyl ring streteh at 244 nm eompared to 785 nm exeitation is shown below. In 
Equation 7.1 all length dimensions are in em, £'Q(244um) + £'^(254nm) is 27500 M'^em'*, 
c is 7.28 M, n (TNT solid) is -1.64, n (TNT solution) is -1.35: 
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When Raman enhaneements for TNT’s phenyl streteh are of the order of mierons 
(10'"^ em), the value of E^ is of the order of . Eor referenee purposes, the deposits 
of explosive residues from fingerprints are typieally less than 7 jum thiek, whieh is 
thieker than the penetration depth of TNT at UV wavelengths [14]. 
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The estimated total Raman enhancement, based on the model presented in 
equation 5.14, at UV wavelengths versus NIR wavelengths for the TNT phenyl ring 
stretch, TNT symmetric NO 2 stretch, and PETN symmetric NO 2 stretch are presented in 
Table 3. 


Raman Band 

Es/D (//m) 

TNT Symmetric NO 2 (229 

750 

nm) 

TNT Phenyl (229 nm) 

2000 

TNT Symmetric NO 2 (229 

740 

nm) 

TNT Phenyl (229 nm) 

6300 

PETN Symmetric NO 2 

70000 

(229 nm) 

PETN Symmetric NO 2 

11000 

(244 nm) 


Table 3. Total Raman enhancements (approximately ±25% error) compared to NIR 

Raman measurements, which factors in frequency enhancements and local field 
effects, are presented for TNT and PETN solids divided by the depth of the 
sample in microns. Eor TNT, the total enhancements for solids are of the same 
order of the resonance enhancements for solutions for micron size samples. Eor 
PETN, which has less laser absorption and Raman signal self-absorption, total 
enhancements are estimated to be greater than enhancements in solution. 

It is important to recognize that this model has not been experimentally verified to 
date and the results from the calculations using equation 7.1 are only estimates. It was 
mentioned earlier that the current model does not include scattering coefficients for the 
solid. In addition, the refractive index values for TNT (solid) and acetonitrile are for 
measurements reported referenced to the Na line (589nm). Refractive index dispersion 
curves for TNT (solid) and acetonitrile are not available in the literature to our 
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knowledge. For TNT, this may not be surprising sinee TNT has several polymorphs that 
would require measurements of each. In general, the refractive index of a material 
increases toward shorter wavelengths, which suggests that the refractive index values 
used in these calculations are on the low side. 

Finally, the local field correction used in these calculations was described in the 
methods section. Briefly, the local field correction has mainly been applied in Raman 
spectroscopy when determining the absolute Raman cross section from liquid 
measurements [37, 38, 45, 47, 48, 49]. The local field correction is added to account for 
intermolecular interactions in the liquid (or solution) phase that are absent in the gas 
phase at low concentration. The local field correction in equation 7.1 is based on that 
used for calculating absolute Raman cross sections for gases using liquid measurements 
and, therefore, may not be the most appropriate correction for accounting for differences 
between solid and liquid phase intermolecular interactions. 
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VIII. CONCLUSION 


A method was demonstrated to estimate pre-resonance and resonance 
enhancements of explosives and select explosive impurities in solutions using acetonitrile 
as an internal standard. Sample transmission and integrated Raman band intensities were 
calculated for sample solutions of PNT; 2,4-DNT; PETN; and 2,4,6-TNT. Spectral 
conditions between NIR Raman and UVRR systems were matched to the greatest degree 
possible; however, the use of the internal standard mitigates differences in spectral 
conditions to allow quantitative enhancement calculations. Enhancement measurements 
were corrected for concentration of the samples, sample transmission, quantum 
efficiencies of the detectors, and detector gain. 

All molecules studied demonstrated appreciable pre-resonance enhancements 
ranging from a low of 20 for PETN’s symmetric NO 2 stretch at 244 nm excitation to a 
high of 2200 for TNT’s phenyl stretch at 244 nm. Additional multiplicative frequency 
enhancements result from UVRR measurements; this includes a 138-fold enhancement 
going from 785 nm to 229 nm excitation and 107-fold enhancement going from 785 nm 
to 244 nm. A model was presented to extend calculated Raman resonance enhancements 
in solution to total Raman enhancements for solid samples. Resonance enhancements are 
observed over a broad bandwidth; the interplay of the optical (i.e., absorption and 
scattering) and physical (e.g., particle size and thickness) properties of the sample and 
sample matrix, the magnitude of the resonance enhancement, and frequency 
enhancements must be explored to find the optimal excitation wavelength. 

Observed polar solvent effects may slightly shift the wavelength where maximum 
resonance effects are observed and also alter the magnitude of the resonance 
enhancement at a given wavelength due to bathochromic and hypsochromic frequency 
shifts. 

In addition to demonstrated Raman signal enhancements, the signal-to-noise ratio 
will also likely be improved in UVRR measurements from the reduction of fluorescence 
interference and other interference due to operation in the UV solar blind region. Einal 
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quantitative assessment of the UVRR standoff Raman signal strengths should also 
ineorporate atmospherie attenuation, espeeially from ozone, and refleetion losses from 
the sample surfaee. 

Sinee different HE moleeules have different exeitation wavelengths where the 
Raman signals are at a maximum intensity, it may be desirable to use a tunable souree for 
explosive deteetion if the required analysis time and speetrometer performanee allows. 
Several deteetion sehemas are possible. A UV wavelength eould be ehosen to give the 
greatest resonanee effeets to a number of prioritized analytes. For example, if TNT, 
PETN and RDX form the most proliferated explosive threats, the eompounds and their 
Raman signal strength at a given wavelength eould be weighted to seleet the best 
performanee of the system. It may also be desirable to seleet two wavelengths that 
optimize performanee for the greatest number of analytes. Another interesting eoneept is 
to eombine UV and NIR or visible measurements to speoifieally look for resonanee 
enhaneements in HE samples to help prevent false positive alarms. 

The overriding goal and motivation of this researeh was to demonstrate possible 
enhaneements that extend to the standoff deteetion of explosives. UVRR speetroseopy 
appears to be more promising in deteeting explosive residues than NIR measurements 
due to resonanee enhaneements, improved seattering effieieney due to frequeney 
enhaneements, and the reduetion in baekground noise. 
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IX. RECOMMENDATIONS 


In order to continue to progress in standoff explosive detection, it is 
recommended that a variety of additional experiments be performed. This research will 
provide additional insight into the military feasibility of Raman spectroscopy in the 
standoff detection of explosives. 

A. HIGH EXPLOSIVE SAMPLES 

TNT, DNT, & PNT were interesting molecules to study to learn more about the 
most important factors that contribute to resonance effects in Raman of high explosives. 
TNT displays stronger resonance enhancements than DNT and PNT. This may be due to 
the combination of increased molecular symmetry and higher molar absorptivity at the 
selected wavelengths of study. The future study of 1,3,5-trinitrobenze, an explosive 
molecule with a high level of symmetry, may help elucidate the factors that are most 
important to resonance Raman. 

Other HE molecules and formulations should also be studied including TATP, 
Composition B, ANFO, C4, Semtex, and nitromethane. 

B. EQUIPMENT DISCUSSION 

The Swept-Wavelength Optical resonance-Raman Device (SWORrD) system 
used in the previous NRL study allowed for the rapid qualitative analysis of Raman 
resonance effects from 210 to 280 nm. A similar tunable system with a 180° collection 
geometry for solid and solution samples would be ideal to rapidly explore Raman 
resonance quantitatively using the method described in this thesis. 

C. MODEL VERIFICATION 

The model presented in Equation 5.14 for estimating the resonance Raman of 
solids using solution data requires experimental verification. In all likelihood, this will 
require basic measurements of some of the unreported optical and physical properties of 
the HE materials. 
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D. STANDOFF DETECTION 


Standoff UVRR measurements using LLNL’s pulsed alexandrite laser system are 
scheduled to begin in the near future under laboratory conditions. The use of a spun 
sample matrix with a solid internal standard should be explored to quantify Raman 
resonance effects at standoff distances. Spinning the sample should help mitigate 
heterogeneity in the sample which can skew results. 

If laboratory standoff UV Raman measurements are favorable, field standoff 
UVRR measurements of a wide variety of target sets would be the next step. 
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APPENDIX A. ADDITIONAL SPECTRA 


A. 785 NM SPECTRA 

1. Standards 


See reference [17] for shift standard peak locations of the Raman shift standards. 



Figure 28. The spectra of acetonitrile, toluene, and a 50% v/v mixture of 
toluene/acetonitrile using a 785 nm laser are presented. 



Figure 29. Raman spectra of cyclohexane used to calibrate the UVRR system at 229 nm. 
Acetonitrile/toluene is not suitable due to the absorption edge in the deep UV. 
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2. Explosives 


The following HE Raman spectra are presented, but were not included in the main 
body of this thesis. 



Raman Shift (cm"^) 


Figure 30. Raman spectra of HMX, ANFO and TATB crystals collected using a 785 nm 
excitation Raman spectrophotometer. ANFO and HMX spectra are offset for 

clarity. 
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Figure 31. Raman spectra of 3,4-DNT; 2,4-DNT; 2,6-DNT; 2,3-DNT; and PNT crystals 
collected using a 244 nm Raman microscope. Laser power was 2.1 mW with 
180- second integration time for all spectra. Spectra are shifted up 2000 counts 
for clarity. Sharp atmospheric O 2 and N 2 Raman lines are visible as indicated. 
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Raman Shift (cm ) 

Figure 32. Normalized Raman speetra of 2,4-DNT crystals collected using at 785 nm (50 
mW), 244 nm (2.2 mW), and 229 nm (0.27 mW) 180 seconds. The symmetric 
NO 2 stretch is the much more intense than the asymmetric NO 2 stretch band, as 
would be expected based on Raman selection rules. The phenyl ring stretch, 
which is very symmetric in nature, is also a strong band. The atmospheric O 2 band 
convolves the NO 2 asymmetric stretch. 
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Figure 33. Raman spectra of 3,4-DNT; 2,4-DNT; 2,6-DNT; 2,3-DNT; and PNT crystals 
collected using a 244 nm Raman microscope. Laser power was 2.1 mW with 180 
second integration time for all spectra. Spectra are offset by 2000 counts for 
clarity. Sharp atmospheric O 2 and N 2 Raman lines are visible as indicated. 



Figure 34. Raman spectra of 9.25 mM PETN in acetonitrile recorded at various laser 
power settings ranging from 20 pW to 2.2 mW using a 244 nm laser excitation. 
Total accumulated CCD integration time for all spectra was 180 seconds to allow 
quantitative analysis of the integrated band intensities to be performed. The 
symmetric NO 2 (1282 cm'*) and (1625 cm'*) and asymmetric NO 2 stretches are 

expanded. 
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Figure 35. Raman spectra of 3.43 mM PNT in acetonitrile recorded at various laser 
power settings ranging from 25 pW to 0.95 mW using 229 nm laser excitation. 
Total accumulated CCD integration time for all spectra was 180 seconds. 



Figure 36. Raman spectra of 2.43 mM 2,4-DNT in acetonitrile recorded at various laser 
power settings ranging from 25 pW to 0.95 mW using 229 nm laser excitation. 
Total accumulated CCD integration time for all spectra was 180 seconds. 
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APPENDIX B. MATLABCODE 


This source code is provided without any warranties on performance. The user is 
advised to test the source code thoroughly before relying on it. The user must assume the 
entire risk of using the source code. 

A. PEAK AREA ALCULATIONS 

% peak_area_v3 

% Created by Christine Paulson, Chance Carter & modified by Maj Bill 
% Short 

% Last modified 03/23/2009 
% Lawrence Livermore National Laboratory 

% This is a MATLAB m-file designed to read a spectra, user select two 
% points for a baseline, compute a baseline, and then integrate the 
% area between the two points. 

close all; clear all; 

cd('C:/Documents and Settings/') % input correct user directory 
% USER SELECTS A FILE 

[FileName,PathName] = uigetfile('*.txtSelect the spectra file.'); 

% LOAD FILE AND CUMULATIVE SUM 
S2 = load([PathName FileName]); 

% Data format: this program will be looking for data in two columns 
% separated by tabs. 

% The X column is comprised of Raman shift data with the y column 
% containing count intensity. 

SI = S2(:,2) '; % select the y column and change the data format to a 
% single row. 

SI = cumsum(SI,1); %cumulative sum, in case of multiple frames 
% MAKE PLOT 

figure; plot (SI, 'k', 'LineWidth', 1); title(FileName); axis tight; 
hold on; 

% USER SELECTS PEAK BOUNDARIES 
datacursormode on 

dispC SELECT THE PEAK BOUNDARIES ON THE FIGURE'); 
dispCHOLD DOWN THE "alt" KEY TO SELECT THE SECOND POINT'); 
quit = 'n'; 
while(quit ~= 'y') 

disp ( ' '); disp ( ' ' ) 

input('Press any key when peak selected.', 's') 
dcm_obj = datacursormode(gcf); 
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info = getCursorInfo(dcm_obj); 

pos = [info.Dataindex]; 

xval = [info.Position] ; 

xval = xval(1:2:end) ; 

disp(['You selected positions: ']) 

disp(xval) 

% COMPUTE BASELINE 
baseline_type = 3; 
xval = sort(xval); 
peakwidth = xval(2) - xval(l); 

if baseline_type == 1 %'midpoint' 

baseline = ones(peakwidth,1).* (SI(xval(1)) + 

SI(xval(2)))/2; 

elseif baseline_type == 2 % 'halfpeak' 

baseline = ones(peakwidth,1).* min([ SI(xval(1)) 

SI(xval(2))]); 

else %baseline linear fit 

increment = (SI(xval(2))-SI(xval(1)))/peakwidth; 
baseline = SI(xval(1)):increment:SI(xval(2)); 

end 

plot(xval(1):xval(2), baseline, 'r', 'LineWidth', 1.5); 

plot(xval(1):xval(2), SI(xval(1):xval(2)), 'r', 'LineWidth', 1.5); 

%COMPUTE AREA 

Sl_subtract = SI(xval(1):xval(2)) - baseline; 
area = sum(Sl_subtract); 

disp(area) 

% Create textbox 

annotation(gcf,'textbox','String',{['Peak Area = ' 
num2str(area)] }, 'FontSize' , 16, . . . 

'FontName','Arial',... 

'FitHeightToText', 'off', . . . 

'LineStyle', 'none ' , . . . 

'Position',[0.1533 0.7667 0.3288 0.1284]); 
quit = input('Do you want to quit? (y/n): ', 's'); 

end 

B. PEAK AREA CALCULATIONS BASED ON HALF OF AN OBSERVABLE 
CONVOLVED PEAK 

% peak_area_v3 

% Created by Christine Paulson & modified by Maj Bill Short 
% Last modified 03/23/2009 
% Lawrence Livermore National Laboratory 

% This is a MATLAB m-file designed to read a spectra, user select two 
% points for a baseline, compute a baseline, and then integrate a peak 
% area of a partially convolved peak by exploiting peak symmetry. 

% User will select the beginning of the peak and peak maximum. Then 
% user selects 'L' or 'R' to tell the 
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% program to compute the full the area from left to right or right to 
% left, based on which portion of the peak is observable. 

%USER SELECTS A FILE 

[FileName,PathName] = uigetfile('*.txtSelect the spectra file.'); 

% LOAD FILE AND CUMULATIVE SUM 
D = load([PathName FileName]); 

%March 20, 2009 - inserted data format change 
SI = D(:,2) ; 

figure; 
plot(D ( :,2)) 

% USER SELECTS PEAK BOUNDARIES 
datacursormode on 

dispC SELECT THE BASELINE POINTS'); 

dispCHOLD DOWN THE "alt" KEY TO SELECT THE SECOND POINT'); 

quit = 'n'; 

while(quit ~= 'y') 

disp ( ' '); disp ( ' ' ) 

input('Press any key when peak selected.', 's') 

dcm_obj = datacursormode(gcf); 

info = getCursorInfo(dcm_obj); 

pos = [info.Dataindex] ; 

xval = [info.Position] ; 

xval = xval(1:2:end) ; 

disp(['You selected positions: ']) 

disp(xval) 

% COMPUTE BASELINE 

baseline_type = 3; 

xval = sort(xval); 

peakwidth = xval(2) - xval(l); 

if baseline_type == 1 %'midpoint' 

baseline = ones(peakwidth,1).* (Si(xval(1)) + Si(xval(2)))/2; 
elseif baseline_type == 2 % 'halfpeak' 

baseline = ones(peakwidth, 1).* min([ Si(xval(1)) 

Si(xval(2))]); 

else %baseline linear fit 

increment = (Si(xval(2))-Si(xval(1)))/peakwidth; 
baseline = (Si(xval(1)) :increment:Si(xval(2))) ; 

end 


cropSl = Si(xval(1) :xval(2 ) ) ; 
figure; 

plot(xval(1):xval(2), baseline, 'r', 'LineWidth', 1.5); hold on; 
plot (SI); 

plot(xval(1):xval(2), cropSl, 'r', 'LineWidth', 1.5); 


figure; 

Sl_subtract = cropSl' - baseline; 
plot(xval(1):xval(2), Sl_subtract, 'g') 
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dispC SELECT THE PEAK POINT'); 
disp ( ' '); disp ( ' ') 

input('Press any key when peak selected.', 's') 

dcm_obj = datacursormode(gcf); 

info = getCursorInfo(dcm_obj); 

pos2 = [info.Dataindex] ; 

xval2 = [info.Position]; 

xval2 = xval2(1:2:end); 

disp(['You selected position: ']) 

disp(xval) 


leftorright = input('Enter "R" for right or "L" 
disp ( ' '); disp ( ' ' ) 
if leftorright == 'R' 

x2 = xval2:xval(2); x2 = x2 - xval(l)+l; 
S1_R = Sl_subtract(x2); 
area_half = sum(Sl_R) 

else 


x2 = xval(1):xval2; x2 = x2 - xval(l) + 
S1_L = Sl_subtract(x2); 

area_half = sum(Sl_L) 

end 


for 


1 ; 


left: 


's'); 


end 


quit 


input('Do you want to quit? (y/n) : 


's' ) ; 
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